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Abstract 

Motivated by the recent observations for the cosmic acceleration and the suitable 
evolution of the Universe provided an interaction (decay of dark energy to matter) 
is incorporated in a cosmological model, we study the cosmological evolution of the 
Interacting Holographic Dark Energy scenario. Critical points are derived and their 
corresponding cosmological models are presented. The dynamical character of these 
models is revealed. 
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1 Introduction 



Recent observations from type la supernovae [1] associated with Large Scale Structure [2] 
and Cosmic Microwave Background anisotropies [3] have provided main evidence for the 
cosmic acceleration. The combined analysis of cosmological observations suggests that 
the universe consists of about 70% dark energy, 30% dust matter (cold dark matter plus 
baryons), and negligible radiation. Although the nature and origin of dark energy are 
unknown, we still can propose some candidates to describe it, namely since we do not 
know where this dark energy comes from, and how to compute it from the first principles, 
we search for phenomenological models. Hopefully, astronomical observations will select 
one of them. The most obvious theoretical candidate of dark energy is the cosmological 
constant A (or vacuum energy) [4, 5] which has the equation of state w = —1. However, as 
is well known, there are two difficulties that arise from the cosmological constant scenario, 
namely the two famous cosmological constant problems — the "fine-tuning" problem and 
the "cosmic coincidence" problem [6]. An alternative proposal for dark energy is the dy- 
namical dark energy scenario. The dynamical dark energy proposal is often realized by 
some scalar field mechanism which suggests that the energy form with negative pressure 
is provided by a scalar field evolving down a proper potential. So far, a large class of 
scalar-field dark energy models have been studied, including quintessence [7], K-essence 
[8], tachyon [9], phantom [10], ghost condensate [11] and quintom [12], and so forth. But 
we should note that the mainstream viewpoint regards the scalar field dark energy mod- 
els as an effective description of an underlying theory of dark energy. In addition, other 
proposals on dark energy include interacting dark energy models [13], braneworld models 
[14], Chaplygin gas models [15], and many others. 

Currently, an interesting attempt for probing the nature of dark energy within the frame- 
work of quantum gravity (and thus compute it from first principles) is the so-called 
"Holographic Dark Energy" (HDE) proposal [16, 17, 18, 19]. It is well known that the 
holographic principle is an important result of the recent researches for exploring the 
quantum gravity (or string theory) [20]. The Holographic Dark Energy model has been 
tested and constrained by various astronomical observations [21, 22, 23, 24] as well as by 
the Anthropic Principle [25]. Furthermore, the Holographic Dark Energy model has been 
extended to include the spatial curvature contribution, i.e. the Holographic Dark Energy 
model in non-flat space [26]. For other extensive studies, see e.g. [27] 1 . 
It is known that the coincidence or, "why now" problem is easily solved in some mod- 
els of holographic dark energy based on this fundamental assumption that matter and 
holographic dark energy do not conserve separately, but the matter energy density decays 
into the holographic energy density [29, 30]. In fact a suitable evolution of the Universe is 
obtained when, in addition to the holographic dark energy, an interaction (decay of dark 
energy to matter) is assumed. 

The remainder of the paper is organized as follows. In Section 2, we present the Holo- 
graphic Dark Energy model where the interaction is included. Since the fractional contri- 
butions of the curvature (Qk), the ordinary matter (f2 m ), and the dark energy (fix) do not 
form a compact state space, we introduce a new set of dimensionless variables which form 
such a compact state space. The evolution equations of these variables are derived and in 
order to study their dynamical character the critical points of this system are obtained. 

1 A very recent development is the idea of bulk holographic dark energy. In this proposal, holographic 
dark energy is accommodated in the framework of braneworld cosmology [28] . 



2 



Cosmological models that correspond to the aforementioned critical points are presented 
and their dynamical character is given. At this point, it is worth noting that one of the 
properties that should characterize a cosmological solution is that the present universe is 
a global attractor, i.e. all the possible initial conditions lead to the observed percentages 
of dark energy and dark matter; once reached, they remain fixed forever [30]. This is the 
reason that here we will focus on the stability of the cosmological models that correspond 
to critical points and more importantly, these cosmological models are attractors. Finally, 
Section 3 is devoted to a brief summary of results and concluding remarks. 



2 Stability of Interacting Holographic Dark Energy 
model solutions 

In this section we consider the Holographic Dark Energy model when there is an inter- 
action between the holographic energy density px and Cold Dark Matter (CDM) which 
has w m = 0. Since we know neither the nature of dark energy nor the nature of dark 
matter, a microphysical interaction model is not available either. However, pressureless 
dark matter in interaction with holographic dark energy is more than just another model 
to describe an accelerated expansion of the universe. It provides a unifying view of differ- 
ent models which are viewed as different realizations of the Interacting Holographic Dark 
Energy Model at the perturbative level [31]. The continuity equations for dark energy 
and CDM are 

Px + 3H(l + w x )p x = -Q, (1) 
p m + 3Hp m = Q . (2) 

The interaction is given by the quantity Q = Tpx with T to have dimensions of pressure. 
This is a decaying of the holographic energy component into CDM with the decay rate T. 
Taking a ratio of two energy densities as r = p m / p x , the above equations lead to 

1+r r 



f = 3Hr 



w x + 



3H 



(3) 



It should be stressed that the decay rate has to be positive, i.e. T > 0, so as to interpret 
a transfer from the dark energy component to the matter component. Moreover, it is 
obvious from equation (3) that any restriction on the parameter of the equation of state 
for the holographic dark energy, i.e. w x , will set constraints on T > [32]. 
Following [33], if we define 

then the continuity equations can be written in their standard form 

p x + 3H(l + w x 1 )p A = 0, (5) 

p m + 3H(l + w c ^)p m = . (6) 
We consider the non-flat Friedmann- Robert son- Walker universe with line element 

1 KT 
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where k denotes the curvature of space with k=0,l,-l for flat, closed and open universe 
respectively. A closed universe with a small positive curvature (f^ ~ 0.01) is compatible 
with observations [34, 35]. We employ the Friedmann equation to relate the curvature of 
the universe to the energy density. The first Friedmann equation is given by 



H 2 + 



k 1 



a 2 3M p 2 



Px + Pm • (8) 



We define the dimensionless fractional contributions of CDM, holographic dark energy, 
and curvature as follows 

m p cr 3M 2 H 2 ' X Pcr 3M*H 2 ' k a 2 H 2 { ' 

Using the results of [36] , the effective parameters of the equation of states given in equation 
(4) are now expressed in terms of the above-mentioned dimensionless parameters as 

wr = _i _ vra; _ (io) 



At this point it should be noted that for the computation of the effective parameters as 
given above, the decay rate is chosen to be 

T = 3b 2 (l + r)H (12) 

where b is a dimensionless coupling constant, and the holographic dark energy density is 
of the form 

p x = 3c 2 M 2 p L~ 2 (13) 

where c is a positive constant which can not be less than 1. The IR cut-off parameter L 
is chosen to be the radius of the event horizon as measured on the sphere of the horizon 
and is defined as 

L = ar(t) (14) 



where r(t) is obtained from 
which yields 



/ 

Jo 



•'■•>_J^ = ^ (15) 

o VI - kr 2 a v ; 

r{t) = -^=wny (16) 

with y = \fkR h /a and is the radius of the event horizon measured in the r direction 

r°° dt 



roo Hf 

R h = a - . (17) 
Jt a 



Now we describe solutions of the model under investigation and also determine their 
stability. To do this, we use a set of convenient phase-space variables similar to those 
introduced in [37, 38]. We introduce the notation similar to those of [37], and work in 
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the 3-dimensional f2-space (Q' k ,Q m ,Qx) where Q! k = —flk- We can now rewrite the first 
Friedmann equation as 



fi A + U, = 1 



It is obvious that when the universe is flat (k = 0) or open (k = — 1), all terms in 
equation (18) that are put together to give unity are nonnegative and take values in 
[0, 1]. Nonetheless all of them are singular in the case of H = 0. For the case of closed 
universe (k = 1) which is of our interest, the terms in equation (18) are not anymore 
all nonnegative. Therefore, the state space defined by the variables (Q' k , Q m , fix) is no 
longer compact (because now Q' k < for k — 1). However, we can employ another set of 
variables describing a compact state space. Therefore, in lieu of using the Hubble function 
H, we introduce the following quantity 



D 



! H 2 



and thus we define the following dimensionless variables 



fix = 



Px 



(19) 

(20) 
(21) 
(22) 



From these definitions it is evident that now the case H = is included. Moreover, the 
Friedmann equation takes the following form 



VL m + Qx — 1 



(23) 



which, together with the fact that — 1 < Z < 1 since the Hubble parameter can be positive 
(H > 0) for an expanding cosmological model or negative (H < 0) for a contracting model, 
implies that the state space defined by the new variables is indeed compact. Introducing 
the primed time derivative 

,_ 1 d 
~ Ddi 

one obtains the system of first-order differential equations 



(24) 



D ' = - z ' D { q+ h) ' 



Z' = z 2 



q + Z 2 [q + — 



where 



—a 



H 2 a 




-3 (l + < ff ) + 2Z 2 [q 
-3 (l + wf) + 2Z 2 (q 

3 (Q x wf + n m w^) 



1 

Z* 
1 



2Z 2 



+ 



2Z 2 



(25) 

(26) 
(27) 
(28) 

(29) 
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The evolution equation for D is not coupled to the rest of the equations (26-28), so we 
will not consider it for our dynamical study. Therefore, we will only study the dynamical 
system for the variables ft = (Z,ft m ,ftx), determined by the equations (26-28). The 
behavior of this system of equations in the neighborhood of its stationary point is deter- 
mined by the corresponding matrix of its linearization. The real parts of its eigenvalues 
tell us whether the corresponding cosmological solution is stable or unstable with respect 
to the homogeneous perturbations [39]. 

To begin with, we have to find the critical points of this dynamical system, which can be 
written in vector form as follows 

ft' = {(ft) (30) 

where / can be extracted from equations (26-28). The critical points, ft*, namely the 
points at which the system will stay provided it is initially at there, are given by the 
condition 

/(fi*)=0. (31) 
Their dynamical character is determined by the eigenvalues of the matrix 



df 
dft 



(32) 



fi=f2* 



The coordinates of the critical points in the state space, i.e. ft = (Z, fl rn , fix), as well as 
their eigenvalues are given in the following table 



Model Coordinates Eigenvalues 

F- (-1, 1,0) -(1 + 3<f, 3<f , 3«f - <)) 

dS_ ("i ' 1 ) ff -(l + 3<,3K ff -< ff ),3<) 

E (°> -sdjF^> 3Kp^y) "(0,-1,1) 

F + (1,1,0) (3<V + 3<,3(<?-<)) 

dS + (1,0,1) (3<,3«-<f),l + 3<) 



with 



a = ^-(l+3<)(l + 3<) . (33 ) 

The cosmological models identified with F + and F_ are the expanding (H > 0) and 
contracting (H < 0) FLRW models (k — A = 0) respectively, dS + and dS- are the 
expanding and contracting de Sitter models (k = p m = 0) respectively, and E stands for 
the Einstein universe (k — 1 and H — 0). 

Concerning the dynamical character of the critical points, if the real part of the eigenvalues 
of a critical point is not zero, the point is said to be hyperbolic [37]. In this case, the 
dynamical character of the critical point is determined by the sign of the real part of the 
eigenvalues: if all of them are positive, the point is said to be a repeller, because arbitrarily 
small deviations from this point will move the system away from this state. If all of them 
are negative the point is called an attractor because if we move the system slightly from 
this point in an arbitrary way, it will return to it. Otherwise, we say the critical point 
is a saddle point. As we can see from the previous table the dynamical character of all 
equilibrium points depend on the equation of state parameters, i.e. and tujf . 



6 



Model 



Repeller 



Attractor 



F_ w5 < -h and wf < wf 



m ^ 3 



< -| and > wf 
£ - 

F+ - wf < -± and wf < wf 

dS+ - wf < -f and wf > wf 

It should be pointed out that for the Einstein universe (E'-model) there is no way to be a 
repeller or, an attractor. It represents a set of infinite saddle points which have to satisfy 
one of the following conditions 

<>~\ and <<~ (34) 
or, wf < ~ and wf > ~ . (35) 

3 Conclusions 

Understanding the dark energy is one of the biggest challenges to the particle physics 
this century. Studying the interaction between the dark energy and ordinary matter will 
open a possibility of detecting the dark energy. Being a dynamical component, the scalar 
field dark energy is expected to interact with the ordinary matters. For example, Carroll 
[40] (see also [41]) has considered an interaction of form QF^F^ with F^ u being the 
electromagnetic field strength tensor which has interesting implication on the rotation of 
the plane of polarization of light coming from distant sources. Recent data on the possible 
variation of the electromagnetic fine structure constant has triggered interests in studies 
related to the interactions between quintessence and the matter fields. 
In this paper we have studied the cosmological dynamics of the Interacting Holographic 
Dark Energy model. We derived the system of the first-order differential equations which 
describes the evolution of the three dimensionless quantities that form a compact state 
space. Furthermore, the five critical points of the aforementioned cosmological model 
were obtained and the dynamical character of these critical points was presented. 
In particular, it was shown that E-model represents a set of infinite saddle points whose 
line element is that of the Einstein universe. By using the fact that the equation of state 
parameters wf ? and should always be negative as indicated by equations (10) and 
(11), we obtained that if 3wf f < 1 the critical points of the F± models are saddle points. 

On the other hand, if 3wf * > 1, the critical point of model F_ will be a repeller under 
the condition |wf | > |wf |, while the critical point of the F + model will be an attractor. 
Concerning the dS± models, we see that if |3w^/ | < 1 then the critical points of dS± are 
saddle points. On the other hand, if |3w^ | > 1, the critical point of model dS_ model 
will be a repeller under the condition |wf ^| < |wf |, while the critical point of dS + will 
be an attractor. 

Finally, it should be stressed that evidence was recently provided by the Abell Cluster 
A586 in support of the interaction between dark energy and dark matter [42]. However, 
despite the fact that numerous works have been performed till now, there are no strong 
observational bounds on the strength of this interaction [43]. This weakness to set strin- 
gent (observational or theoretical) constraints on the strength of the coupling between 
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dark energy and dark matter stems from our unawareness of the nature and origin of 
dark components of the Universe. It is therefore more than obvious that further work is 
needed to this direction. 
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